Preprint ICN-UNAM 16/04 



Invited contribution to Czechoslovak Chemical Communications 
to a Special Issue in honor of Professor Josef Paldus 

Hydrogen atom and one-electron molecular systems in a strong 
magnetic field: are all of them alike? 

A. V. Turbinery A. B. Kaidalovjj] and J. C. Lopez Vieyrjf] 
Instituto de Ciencias Nucleares, Universidad Nacional Autdnoma de Mexico, 
Apartado Postal 70-543, 04510 Mexico, D.F., Mexico 
(Dated: February 2, 2008) 

Abstract 
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The behavior of Coulomb systems in a strong magnetic field has always attracted a lot 
of attention. This was justified b y th e presence of strong magnetic fields in astrophysics 
(neutron stars and white dwarfs 26]), as well as in plasma and semiconductor physics. 
In particular, for many years it has existed a question about the content of the neutron 
star atmosphere. Since the seminal papers by Kadomtsev-Kudriavtsev |l| and Ruderman 
[2I it was believed that the neutron star atmosphere subject to a strong magnetic field is 
made from atomic-molecular compounds. However, in order to construct a model of the 
atmosphere it is necessary to explore matter and its properties in a strong magnetic field. 

Recently, it was discovered that the interplay of Coulomb and magnetic forces for B > 
10 11 G leads to a new physics: new bound one-electron Coulomb systems appear like the 
exotic molecular ions fhj~ + [3], (HeH) ++ and He 3 ^ pf. These Coulomb systems do not 
exist without magnetic field. All of them are characterized by very large binding energies 
growing with magnetic field. For all magnetic fields, where the non-relativistic considerations 
are justified (B < 4.414 x 10 13 G), of one-electron atomic-molecular systems the hydrogen 
atom - the only neutral system- is characterized by the highest (!) total energy, being 
correspondingly the least bound system At the same time it seems natural to assume 
that more-than-one-electron Coulomb systems in a strong magnetic field are not strongly 
bound (if bound) due to the fact that all electron spins should be parallel, being antiparallel 
to the magnetic field direction. 

All one-electron molecular systems have a certain common feature. Their optimal con- 
figuration is always a configuration where all massive charged centers are situated on a 
magnetic line. We call it parallel configuration. Only these configurations are considered in 
the present article. 

It is well known that studies in a strong magnetic field are very complicated for several 
reasons. Perhaps, the most serious, conceptual reason is related to the fact that the bound 
states are of a weakly-bound-state nature (the binding energies are much smaller than the 
total ones). The perturbation theory in powers of B is fast divergent and thus cannot be 
used. Asymptotic expansions at B = 00 have extremely complicated form but, usually, 
have no domain of applicability inside non-relativistic considerations. For example, it can 
be easily checked that at an extremely strong magnetic field near the edge of applicability of 

n 

non-relativistic approximation B = 10000 a. u.(= 2.35 x 10 13 G) [27], the ionization energy 
Eb calculated numerically differs by 300% (!) from the value obtained using the leading 
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^asymp = 



term in asymptotic expansion 



where is given in atomic units as well as B (see [6]). Another example is given by the 
ratio of the binding energies of H^" and Hjj~ + . Asymptotically, for B tending to infinity, 
this ratio should be 1/2.25 . However, the numerical results at B = 3 x 10 13 G p] give 
« 1 for this ratio. Therefore the only methods which can be used are either numerical 
or variational. For B > 10 11 G, to the best our knowledge, the numerical methods were 



(see , for example, the excellent early review the book 



used for the hydrogen atom on 

and recent review articles 9|, llOj). Usually, these methods are very slow-convergent 
and extremely difficult to implement. The most popular method to study one-electron 
molecular systems is the variational method (for review, see (Uj). However, the use of the 
variational method is associated with a difficult procedure of minimization and, sometimes, 
with numerical calculation of multidimensional integrals with high accuracy, which can also 
be quite cumbersome. In any case, the calculations are made for some particular values 
of magnetic field. It seems natural to create some approximate expressions valid for all 
magnetic fields, even having not high accuracy, in order to make at least rough estimates. 

The accurate results of calculations of different quantities for low-lying states reveal a 
smooth, simple-looking behavior with rather slow changes with magnetic field. However, a 
straightforward attempt to construct approximations either fails or leads to quite compli- 



cated expressions, at least, at first sight (see, e.g., |l2j,|l3j). Physical intuition gives a feeling 
that there must exist a certain qualitative technique, for example a type of semi-classical 
approximation providing an approximate qualitative description of these results. So far it 
is not clear how such a technique can be approached. A goal of this paper is to consider a 
certain simple alternative to this unclear-how-to-approach technique - to build approxima- 
tions of the main characteristics of the one-electron atomic-molecular systems in a constant 
uniform magnetic field in their lowest state, such as total and binding energies, equilibrium 
distances, electron cloud sizes, quadrupole moment, by following simple physical arguments. 
Our basic assumption is that the physics is mainly governed by a single parameter: the ratio 
of the transverse to longitudinal size of the electron cloud. Of course, in dimensionful quan- 
tities such as equilibrium distances or quadrupole moment the transverse and longitudinal 
sizes should appear explicitly but only in a form of parameters which carry a dimension. 
Hereafter we denote the transverse size of the electron cloud as r t , and the longitudinal size 



as r\. 

As always we consider the one-electron Coulomb systems with infinitely-heavy charged 
centers, protons and/or a-particles (the Born-Oppenheimer approximation of the zero order) 



situated on the z-axis 



281 ] . If these charged centers are of the same charge, they are assumed 



to be identical. Although we use the word 'proton' it implies that in the Born-Oppenheimer 
approximation it can be deuteron or triton. The magnetic field of strength B is directed 
along the z axis, B = (0,0, B). Throughout the paper rydberg (Ry) is used as the energy 
unit. For the magnetic field we use either atomic units or Gauss (G) with the conversion 
factor Bq = la.u. = 2.35 x 10 9 G. For the other quantities standard atomic units are used. 
The distances between infinitely-heavy charged centers are denoted by R letters, whereas the 
distances between centers and electron are denoted by r letters. The distance between the 
electron position and the z axis is denoted by p. In particular, the potential corresponding 
to the hydrogen atom is given by 



where p — y x 2 + y 2 and r is the distance from the electron to the charged center. The 
potential 

_ 2Z X Z 2 _ 2Z X _ 2Z 2 B 2 p 2 
V ~ R n r 2 + 4 ' {Z) 

describes the ions (the system (ppe), Z\ = Z 2 = 1), (HeH) ++ (the system (ape), Z\ = 
1,Z 2 = 2), He2 3+ ' ) (the system (aae),Zi = Z 2 = 2), where r\{r 2 ) is the distance from the 
electron to the charged center 1 (2) and R is the distance between the charged centers. In 
turn, the system Hjj~ + is described by the potential 

2 2 2 2 2 2 B 2 p 2 

~R~l + ~r1 + r + + rZ~ n~ 7 2 ~ 7 3 + 4 ' ^ 

where Ti is the distance from the electron to the charged center i and R± are the distances 
from the central charge, placed in the origin, and the side charged centers. The equilibrium 
distance, which corresponds to the minimum of the total energy, is defined by the distance 
between the most-distant charged centers, which is L eq = R eq for the two-center case of 
H^", (HeH) ++ , He2 + and L eq = R +eq + R_ cq for three-center case of Hg~ + . Hereafter, magnetic 
field is defined in dimensionless units (a.u.) as B/Bq, where Bo = 2.35 x 10 9 G, which we 
continue to denote as B. 
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I. THE H-ATOM 



Let us take hydrogen atom - the simplest one-electron system - placed in a constant 
uniform magnetic field B directed along the z-axis. Due to the Lorentz force, the spherical 
symmetric electron cloud (in the absence of a magnetic field) is deformed to a cigar-like 
form. The size of the electron cloud r t in transversal direction to z-axis shrinks drastically 
~ B~ l l 2 at large magnetic fields, being close the value of the Larmor radius. As to the 
longitudinal size r\ it also contracts at large magnetic fields but at a much more moderate 
rate, ~ (log-B)^ 1 (see, e.g., J2, 14, An interplay of these two types of behavior explains 
the cigar-type form of the electron cloud. At very large magnetic fields, the cigar-type 
form evolves to a needle-like form known as the Ruderman needle. In Fig. la the form 
of the electron cloud is illustrated for B = 10 12 G j^]. In particular, the longitudinal 
size of the electron cloud shrinks in comparison with the zero-magnetic-field case about 
four times. Therefore, the apparent classical (electrostatic) appearance of the magnetic 
field influence is characterized by a change of the form of the electron cloud, which can be 
roughly approximated by the ratio of two classical parameters r t ,r\. In fact, it is the major 
assumption of the present approximation scheme. We also assume that these parameters 
r t , ri are defined by the expectation values, 

r t =< p > , n = 2 < \z\ > . (4) 

If a definition of the transversal size r t looks natural from the physical point of view and 
rather unambiguous, definition Q of the longitudinal size is not so obvious. It can be chosen 
as ~ 2 2 >, or as a linear combination of < \z\ > and V< z 2 >. So far it is not so clear 
what would be physical arguments which allow to specify a definition. Eventually, it turns 
out it is not very important what quantity is used to define r\. The results of the fit remain 
very similar although there can be some difference in the parameters. 

The binding energy is by definition the difference between the energy of free electron 
in magnetic field (the Larmor energy) B and the total energy of the atom, = B — E^. It 
is known that E^ in the weak-field regime is represented by the Taylor expansion in powers 



n 



of B 2 , while for large B it behaves ~ (log 2 B) (see, for example, [bjj and discussion in Q|). 
Following the above assumption the binding energy depends on the ratio X = r t /r\, 

E h = E h (X) . (5) 
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FIG. 1: The contours of the electron clouds on the plane x = for the H atom (a), ion (b) 
and H^" + (c) placed in a magnetic field B = 10 12 G directed along the z-axis, (x, y, z) are in a.u. 
For the H^" , H^-ions the protons are situated on the z-line. The values of 2 < \z\ > are shown 
by bars: the short bars correspond to B = 10 12 G , they are compared with the long bars which 
correspond to B = for H, and to B = 10 10 G (near the threshold of existence) for H^ + . It 
illustrates a shrinking of electronic longitudinal size with magnetic field growth. 



It is quite natural to approximate the transverse size r t =< p > as follows 

r^^4^rJ\^), (6) 



(l + a t 2 £ 2 )V4 \l + b t B 2 / 

where r^,a t ,a t ,b t are parameters, which are found by fitting the calculated expectation 
values for < p >. The formula (J0J) is written in such a way as to reproduce a functionally- 
correct perturbative expansion of < p > at B = (in powers B 2 ) and r® = v^ob? where 
a-Q — 1 a.u. is the Bohr radius. At large B the right-hand side of Eq.© behaves as ~ B~ x l 2 
simulating the Larmor radius behavior. 

In Fig.Elone can see that Eq.© fits data on < p > calculated using the formalism developed 
in 12j with accuracy better than one percent at B > 10 9 G. The parameters of the fit are 
given in Table HJ The parameter is also found from the fit, it deviates from \[2 (see 
above) by ~ 8%. It reflects the fact that the accuracy provided by formula (JBJ) diminishes 
as the magnetic field decreases (see the discussion below). 

At first sight, it is a much more complicated task to describe the longitudinal size, 
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FIG. 2: H atom: the fit of transverse size of the electron cloud < p > versus B using Eq.©. 
Calculated values are indicated by crosses. 



TABLE I: The parameters of the fit (jHJ) of the transversal size of the electron clouds of the H atom 
and H^,H^ + ,(HeH)++,He^ + ions (in a.u.). 



System 


1 1 


«t 


at 


h 


H atom 


1.17533 


0.44904 


1.20981 


1.81098 


H+ 


0.954427 


0.23615 


0.376237 


0.62194 


H 3 ++ 


0.645875 


0.048196 


0.00970609 


0.0230488 


(HeH)++ 


0.174416 


0.019657 


0.00000030 


0.0000003 


He|+ 


0.200825 


0.026449 


0.00000150 


0.00000148 



r\ = 2 < \z\ >. The approximation we propose to use is 



n 



1 + aiB^ 



(7) 



1 + oi log(l + ffB 2 + 7l 2 S 4 ) V 1 + b y B 2 

where r®, at\, fti, 71, a\, b\ are parameters, which are found by fitting the calculated expectation 
values for 2 < \z\ >. Formula (0) has the perturbative expansion in powers B 2 , which agrees 
with perturbation theory results and rf = 3/2 qb, where = 1 a.u. is the Bohr radius. At 
large B, the right-hand side (J7J) behaves as ~ (log B) 1 as should be in accordance with the 
qualitative arguments. 

n Fig. Elone can see that (JJJ) fits data on 2 < \z\ > obtained in the formalism developed 
in 3| with accuracy better than 1% at B ~> 10 9 G. The parameters of the fit are given in 
Table HU The parameter rf is also found from fit. Surprisingly, it deviates from 3/2 (see 
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FIG. 3: H atom: the fit of longitudinal size of the electronic cloud 2 < \z\ > versus B using Eq.(|7|). 
Calculated values are indicated by crosses. 
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FIG. 4: H atom: ratio X = 2 <f^|'> • Calculated values are indicated by crosses. 

above) insignificantly, by < 1 %, in contrast to what happened for the parameter r°. 

In Fig. |3] a comparison of the ratio X = r t /r\ (see Eqs. ©-©) with parameters taken 
from Table H] with results of calculations is presented. One can clearly see that both data and 
fitted curves of Figs. 01 and ^demonstrate a certain irregularity in the range (5 — 50) x 10 10 G. 
It is a transition region from the Coulomb regime, where the Coulombic forces dominate over 
magnetic forces to the Landau regime where in the (x, y) plane Coulombic forces become 
sub dominant. 

Following the assumption (JHJ let us approximate the binding energy 

E b = AX^ + BX Y + C , X x = \ogX , (8) 
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TABLE II: The parameters of the fit of the longitudinal size of the electron clouds of H atom, 
IFj~,IF|~ + molecular ions as well as (HeH) ++ ,He2 + using by (JTJ) (in a.u.). 



System 


~0 

r i 


a\ 


A 


71 




h 


H atom 


1.49719 


0.179332 


0.320252 


0.001164 


1.07512 


1.35162 


H+ 


1.72041 


0.254255 


0.141004 


0.0004436 


0.340131 


0.497807 


H 3 ++ 


1.94408 


0.279956 


0.0191558 


0.000008 


0.0066712 


0.0136894 


(HeH)++ 


1.72219 


1.155934 


0.405312 


0.000311 


0.131457 


0.0600481 


Hel+ 


0.65727 


0.228616 


0.0011334 


0. 


0.0000116 


0.0000228 




B (a.u.) 

FIG. 5: H atom: the fit of the binding energy using (|8*j). Calculated values are indicated by crosses. 



where A, B, C are parameters, which are found by making fit of the results of calculations 
of the binding energy. It is worth emphasizing that the parameters of X(X{) are already 
fixed by the fits (|B |L (J7 j) of r t , ri, respectively. The formula (JHJ) agrees with the perturbative 
expansion in powers B 2 (at small B) and gives a correct asymptotic expansion at large B. 



i of the best known results for the 
121 ]. The parameters A, B, C are given 



In Fig. El it is shown the fit using the formula 
binding energies from combined with those from 
in Table ILTT1 In the whole range of explored magnetic fields 10 9 — 4.414 x 10 13 G, formula 
(JHJ) approximates the binding energies with a relative accuracy, which does not exceed few 
percent and becomes more accurate with growing magnetic field. It is worth mentioning 
that when the parameter A = 4 in the approximation (JHJ) its asymptotic coincides with the 
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exact asymptotic (see, e.g., [14|, §112) 



E h « log 



B 



B 



oo 



(9) 



where is in Ry. In fact, the deviation [A/A — 1| gives a feeling about the quality of our 
approximation. Clearly, this estimate is very rough ca. 20% (see Table III), while a real 
accuracy of approximating the binding energy is a few percent. 

TABLE III: The parameters of the fit of the binding energy of the H atom, H^,H^ + molecular 
ions as well as (HeH) ++ ,He2 + using (jSJ) (all in Ry). 



System 


A 


B 


C 


H-atom 


3.22532 


0.53945 


1.37932 


H+ 


8.23442 


6.8246 


2.99945 




12.8455 


20.4849 


3.95821 


(HeH)++ 


15.7401 


6.1134 


-5.3756 


Hel+ 


26.2926 


32.9181 


-0.28129 



One of the important characteristics of the magnetic field influence on the H-atom is the 
appearance of the quadrupole moment 



Q = -Q z 



2 < z 2 > - < p 2 > . 



(10) 



Recently, the first quantitative study of the quadrupole moment was carried out 
formula (fTUj) suggests immediately the following approximation 

Q = 2r 2 (A q - a 9 Xi) - r 2 (B q + b q X x ) , X x = log X , 



121. The 



(11) 



where A q = 0.325447, a q = 0.049432,5, = 1.32012,6, = 0.955362 are dimensionless param- 
eters, which are found by fitting the quadrupole moment. The parameters of X\ are already 
fixed in the fit of parameters r t , r\ using © and (J3), resp ectively. Formula (jllj) describes 
correctly the expansion at small and large B (see [2], [l2j, [l7J ) . It fits the results of calculations 
in [l^ with an accuracy of few percents (see Fig. (HJ) - 

We made an analysis of the expectation values < |z| n > at n = 2,3,4,5. It turns out 
that the calculated expectation values admit a very accurate polynomial approximation in 
terms of a single expectation value < \z\ >, 



<\z\ n >= P n (<\z\ >) 
10 



(12) 
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FIG. 6: H atom: the fit of the quadrupole moment using (fTTj) . Calculated values (see Q]) 
indicated by crosses. 



are 



where P n is a n-th degree polynomial. It seems natural to assume that ()12|) holds for any 
n, hence any expectation value is defined by < \z\ >. This leads to a striking hypothesis 
that the ground state eigenfunction integrated over p can be viewed as a one-parametric 
probability distribution (!). 



II. THE H7, MOLECULAR ION 

In this Section we consider the molecular ion in parallel configuration, when the 
protons are situated along the magnetic line. The form of the electron cloud is illustrated 
in Fig. lb for the magnetic field B = 10 12 G. The transversal size of the electron cloud r t 
shrinks drastically, ~ B^ 1 ^ 2 , at large magnetic fields, being close to the value of the Larmor 
radius similarly to what happens for the hydrogen atom. As to the longitudinal size n it 
also shrinks but at a much slower rate ~ (logi?) -1 . In particular, the longitudinal size of 
the electron cloud shrinks in comparison with vanishing magnetic field about five times (see 
Fig. lb). 

Following the same arguments which were used earlier for the H atom, we again assume 
that the dynamic characteristics of the in a magnetic field depend on the expectation 
values of transversal (r t ) and longitudinal (n) sizes. The dependence of them on magnetic 
field is approximated by similar formulas The binding energy at equilibrium distance 

between protons depends on the ratio X = r t /r\. Eventually, the binding energy is written 
in the same form (jHJ) with the same expressions © and (J7J) as is done for H atom but with 
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FIG. 7: ETj~ ion: a fit of the transverse size of the electron cloud < p > using Eq.©. Calculated 
values are indicated by crosses. 
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FIG. 8: ion: a fit of the longitudinal size of the electron cloud 2 < \z\ > using Eq.Q (solid 
curve) and of the equilibrium distance R cq (dashed curve) using Ea. ((T^)) . Calculated values are 
indicated by crosses. 



different parameters. For the fit we use the results of recent calculations of the binding 



energy which were carried out in [19|, |2Cj. These parameters of the fit are presented in 
Tables HMTTTI and the fit is illustrated by Figs. CHI 

In order to approximate the equilibrium distance R eq we assume that R eq is proportional 
to the longitudinal distance n with a small correction in X\ 

r° ( \ + a\B 2 \ 

ft, = r,(c„ + + ctf) = TT -^ Tm (^ j (co + c lXl + ctf) , (13) 

where the parameters Cq,c\,C2 are found from the fit of the results of calculations of the 
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FIG. 9: ion: the ratio X = 2 <\^\> • Calculated values are indicated by crosses. 




B (a.u.) 

FIG. 10: H,j~ ion: the fit of the binding energy using (|SJ). Calculated values are indicated by crosses. 

equilibrium distance which were carried out in jigL |^ . The parameters of the fit are given 
in Table IIVI The fit is illustrated in Fig. |HJ It is worth mentioning that the parameters 
Cj, i — 0, 1, 2 decrease very fast with i (see Table HVj) . This can be considered as an indication 
of adequateness of the approximation formula (|13|). 

Similarly to what happened for H atom, the plot of the ratio X = r t /r\ (Fig. |UJ) reveals 
a certain irregularity in behavior of the calculation results as well as the fit in the range 
(5 — 50) x 10 10 G. We assign these irregularities to a transition from the Coulomb to the 
Landau regime. An overall quality of the fit for the domain 10 9 — 4 x 10 13 G is very high, 
about 1-2% except for the above-mentioned region where the accuracy drops to 5-10%. 

Similar to what was done for the H-atom we carried out a calculation of expectation 
values < \z\ n > , n = 2,3,4,5. It turns out that these expectation values admit very 
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TABLE IV: The dimensionless parameters of the fit ()13)) of the equilibrium distance (in a.u.) of 
H+,H++ and (HeH)++, Herons . 



System 


CO 


Cl 


C2 


H+ 


1.37384 


0.389879 


0.0430844 


H 3 ++ 


4.48200 


2.25814 


0.380948 


(HeH)++ 


4.15754 


2.31113 


0.409048 


Her 


1.83774 


0.51165 


0.0626179 



accurate polynomial approximation in terms of the expectation value < \z\ > (see (112)) ). It 
seems natural to assume that ()12|) holds for any n. This leads to the hypothesis that the 
ground state eigenfunction integrated over p defines a one-parametric distribution similar to 
what appears for the H atom (see previous Section). 



III. THE H++ MOLECULAR ION 

Now we consider the exotic system hht + theoretically predicted in |3| , which is made out 
of three protons situated along the magnetic line and one electron (parallel configuration). 
This system appears as a quasi- stationary state at B > 10 10 G J^. The form of the electron 
cloud for B = 10 12 G is shown in Fig. lc. It is clearly seen that the transversal size of the 
electron cloud r t shrinks drastically, ~ B~ l l 2 at large magnetic fields similar to what happens 
for the hydrogen atom and the molecular ion which is of the order of the Larmor radius. 
As to the longitudinal size r\ it also contracts but in much slower rate, ~ (log-B)" 1 , at large 
magnetic fields. 

We follow the same idea of approximation as for H and assuming that the physics 
is governed by a single parameter X = r t /r\. The same approximation formulas © and 
(J7J) are used for the transverse (r t ) and longitudinal (r\) sizes, respectively, as it is done for 
H-atom and H^. Their parameters are found by fitting the results of calculations. The data 
for r t ,ri are obtained using a strategy described in 2fJ. The parameters of the fit are given 
in Tables GJ |TI] The fit of r t and r\ is illustrated in Figs. [TO [TO Fig. demonstrates the 
behavior of the X . The binding energy Eb which is calculated in |2Cj is approximated using 
the formula (jSJ) (see Table llTH for parameters of the fit). The fit is illustrated in Fig. El 
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B (a.u.) 

FIG. 11: ion: the fit of transverse size of the electron cloud < p > using Eq.ljB]). Calculated 

values are indicated by crosses. 

In the same way as it is done for H^, we assume that the equilibrium distance between 
protons are mostly defined by the longitudinal size of the electron cloud (see (J7J)), which are 
slightly modified by including the terms depending on X\ = logA\ Finally, the equilibrium 
distance is approximated by Eq. (jl3|) as was done for (see Fig. I12p. The parameters of 
the fit are given in Table ITVl It is worth mentioning that the parameters % = 0,1,2 
decrease very fast with i. This might be considered as an indication of adequateness of the 
approximation formula (|13jl . 

In the fit, some irregularities can be seen in the region (5 — 50) x 10 10 G, near the threshold 
of appearance of the FFj~ + ion (see Figs. [TT1 - IT3jl similarly to those that were observed for 
the H-atom and for the H^~-ion. One of the reasons for these irregularities can be related to 
highly increased technical difficulties we encountered exploring this region. This could lead 
to a loss of accuracy. The overall quality of the fit for the range 10 11 — 4.414 x 10 13 G is very 
high, 1-5%. 

Similarly to what was done for the H atom and the H^" molecular ion, we calculate the 
expectation values < \z\ n > , n = 2,3,4,5 for the ion. It turns out that these expec- 
tation values admit a very accurate polynomial approximation in terms of the expectation 
value < \z\ >, see Eq. (112)1 . It seems natural to assume that (|T2*jl holds for any n. The ground 
state eigenfunction integrated over p seems to define a certain one-parametric distribution. 
A similar phenomenon occurs for the H atom and the molecular ion. 
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FIG. 12: Hjj~ + ion: the fit of longitudinal size of the electron cloud 2 < \z\ > using Eq.((7|) (solid 
line) and of the equilibrium distance L eq = 2i? eq using Eq, (|l^|) (dashed line). Calculated values 
are indicated by crosses. 




IV. THE (HeH) ++ MOLECULAR ION 

Recently, it was theoretically predicted that the exotic molecular ion (HeH) ++ can exist 
for B > 10 12 G Q. Following the same idea of approximation as it was implemented for the 
H atom and for the H^", H^~ + molecular ions, we can construct high-accuracy approximations 
for the exotic (HeH) ++ ion. Transversal (r t ) and longitudinal (n) sizes [30] of the electron 
cloud as a function of the magnetic field are approximated by the expressions (JHJ) and (JJJ) 
(see Figs. [El and Uni). The parameters of the approximations ©-© obtained through fitting 
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FIG. 14: Hj|~ + ion: the fit of the binding energy using (|5|). Only calculated values which are 
indicated by crosses are used for fitting, while calculated values shown by circles are not taken into 
account (see text). 
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FIG. 15: (HeH) ++ ion: the fit of the transverse size of the electron cloud < p > using Eq.(jHJ). 
Calculated values are indicated by crosses. 



the data from |4| are presented in Tables H1HT1 In Fig. El the ratio X is compared with the 
calculated data from |4|. The fit of the binding energy was performed using the formula 
(JBJ) (see Fig. [TSj) . The parameters of the fit are presented in Table ITTT1 For the equilibrium 
distance R eq , the approximation (fl"3|) is used (see Fig. IT6|) : the parameters are presented in 
Table |TV1 The overall quality of the fit for the range 10 12 — 4.414 x 10 13 G is very high, 
around 1 %. 



17 




1000 10000 
B (a.u.) 



FIG. 16: (HeH) ++ ion: the fit of longitudinal size of the electron cloud < (z — z ma _ x ) > using Eq.Q 
(solid line) and the equilibrium distance R eq using Eq.([13[) (dashed line). Calculated values are 
indicated by crosses. All data in a.u. 
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FIG. 17: (HeH) ++ ion: the ratio X = <p> w . Calculated values are indicated by crosses. 

V ' <(2-Zmax)> J 



V. He| + MOLECULAR ION 

Recently, it was theoretically predicted that for B > 100 a.u. the exotic molecular 
ion He2 + can exist j^. Following the same idea of approximation as for the H atom and 
the H^, H^" + , (HeH) ++ molecular ions (see previous Sections), we would like to construct 
accurate approximations for the exotic He^ ion. Transversal (r t ) and longitudinal (n) sizes 
of the electron cloud as a function of the magnetic field are approximated by the expressions 
O and O, respectively (see Fig, M- M- The parameters of the approximation, © 
and (0) obtained through the fit of the data obtained in |4( are presented in Tables HI- HTJ 
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FIG. 18: (HeH) ++ ion: the fit of the binding energy using Calculated values are indicated by 
crosses. 



respectively. In Fig. the ratio X is compared with the calculated data from |4j . The fit 
of the binding energy was performed using the formula (JHJ) (see Fig. I22J1. The parameters of 
the fit are presented in Table ITTT1 For the equilibrium distance R eq the approximation (JT3J) 
is used (see Fig. EOJ) with parameters presented in Table II VI 

Some irregularities can be seen in the fit in the region (2 — 5) x 10 11 G, near the thresh- 
old of appearance of the He2 + ion (see Figs. ITHl l2*2*|) similar to those which were observed 
for the H atom and for the H^~,lHLt + ions. The overall quality of the fit for the region 
10 12 - 4.414 x 10 13 G is very high, around 1 %. 
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B (a.u.) 

FIG. 19: He 3 ^ ion: the fit of the transverse size of the electron cloud < p > using Eq.JHJ). Only 
calculated values indicated by crosses, are used for fitting, while the calculated values shown as 
black squares are not taken into account (see text). 



VI. CONCLUSION 



In this work we presented a phenomenological model of the behavior of different one- 
electron atomic-molecular systems in a strong magnetic field. The model is based on a 
surprisingly simple physical idea that the ground state depends on a ratio of transverse 
to longitudinal size of a system placed in a strong magnetic field only. Since accurate 
numerical studies in a strong magnetic field are very tedious from a technical point of view a 
construction of a phenomenological model which provides approximate expressions for basic 
characteristics of a system for any value of a magnetic field strength can be quite useful for 
applications. 

One of the motivations of the present work is related to the fact that the neutron star 
atmosphere is characterized by strong magnetic fields, 10 12 — 10 13 G. It seems natural to 
anticipate a wealth of new physical phenomena there. However, for many years the ob- 
servational data did not indicate anything unusual, corresponding to the black-body radi- 
ation. On 2002, the CHANDRA X-ray observatory collected data on an isolated neutron 
star 1E1207. 4-5209 which led to the discovery of two clearly-seen absorption features at 
~ 0.7 keV and ~ 1.4 keV [21]. It is necessary to mention that the XMM-Newton X-ray 
observatory recently confirmed the results of Chandra/ ACIS related to absorption features 
at 0.7 and 1.4 keV |22|. We proposed a model of hydrogen atmosphere with main abundance 
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B (a.u.) 

FIG. 20: He^" 1 " ion: the fit of longitudinal size of the electronic cloud 2 < \z\ > using Eq.(|7jl 
(solid line) and of the equilibrium distance R eq using Ea. (|13|) (dashed line). All data in a.u. Only 
calculated values which are indicated by crosses are used to make a fit, while calculated values 
shown by triangles (for longitudinal size data) and by black squares (equilibrium distance data) 
are not taken into account (see text). 
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FIG. 21: He^ ion: the ratio X 



2<|z|>- 



of the exotic molecular ion which explains these absorption features assuming that the 
surface magnetic field is ~ 5 x 10 14 G j^j]. For other neutron stars, observational indications 
of the existence of absorption lines in their spectra were already found Q, Q- 

It seems 

natural to anticipate forthcoming observations of other neutron stars which will likely reveal 
absorption features. The study presented here can be of certain use in identifying possible 
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FIG. 22: He^ ion: the fit of the binding energy using (jBJ). 

absorption features. 
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